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ABSTRACT 



Three numerical model simulations of cyclogenesis are compared to examine the role 
of boundary layer stratification in enhancing the components forcing the ageostrophic 
circulation of an idealized warm front. The Sawyer-Eliassen diagnostic equation is ap- 
plied to examine the contributions of frictional forcing and diabatic heating, as well as 
confluent and shear geostrophic deformation, to forcing of the sccondaip- circulation of 
a warm front. The surface heat and moisture flux distributions are varied in each case, 
in order to evaluate the effect on each component as well as the cyclogenesis. 

Results confirm previous studies that geostrophic deformation forces strong 
frontogenesis at the surface, and at mid-levels frontogenesis is weaker and forced pri- 
marily by latent heat release. Trontogenetical forcing is modified by small-scale fric- 
tional forcing and diabatic heating, which depend upon the surface and boundarx’ layer 
processes. Although frictional forcing comprises less than five percent of total forcing 
of the warm front, it apparently enhances Irontogenesis. partially due to indirect eflccts 
on other more dominant processes. The intensity of frictional forcing is strongly de- 
pendent on the surface heat flux distribution and track of the cyclone relative to the 
sca-stirface temperature gradient. In the absence of surface fluxes, frictional forcing is 
negligible. Surface forcing due to diabatic heating is frontolytical. and reduces total 
forcing at low levels by about twenty percent. The magnitude of the frontolytical forcing 
by diabatic heating at low levels is only partially dependent on the surface heat and 
moisture flux distribution. Even in the absence of surface heat and moisture fluxes, the 
frontolytical forcing persisted, suggesting problems with the cumulus parameterization 
or errors in the computation of this term. These results remain to be verified against an 
actual case study of a warm front. 
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I. INTRODUCTION 



A. CVCLOGENESIS AND FRONTOGENESIS 

Mid-latitude cyclones are transient waves in the global zonal flow that perform an 
important role in transporting heat and momentum in the atmosphere. Historically, 
cyclone development was explained as a wave instability on the polar front, with the 
existence of the front considered a necessary precursor for cyclogenesis. This theoiy, 
known as the polar front theoiy of cyclogenesis, significantly influenced Petterssen 
(1956) who related cyclogenesis to the interaction of vorticity advection aloft with sur- 
face frontal features. The modern view, based on baroclinic instability theoiy, considers 
fronfj to be a consequence of the development of a cyclone and a source of baroclinic 
energy for cyclogenesis. The baroclinicity of the atmosphere and the synoptic-scale flow 
combine to lead to cyclogenesis, the development of sub-synoptic scale fronts, and the 
lormation of jet streaks. 

1 he basic premise of baroclinic instability theoiy is that the circulation of a cyclone 
1 C .ills from any small perturbation forced by advection of positive (cyclonic) vorticity 
aloft with the advance of the upper-level wa^ e or by warm ad^'ection at low levels over 
a broad region. Geostrophic deformation associated with cyclogenesis acts to enhance 
the temperature gradient in particular regions of a cyclone. Ihis process is known as 
Ironiogcnesis. Vertical and divergent ageostrophic motions arise in response to 
Irontogenesis. and are ihe mechanism for the adjustment of vorticity and thermal prop- 
erties to return the system to thermal wind balance. The ageostrophic circulations as- 
sociated with frontogencsis are small-scale manifestations of the basic synoptic-scale 
forcing. .Although these synoptic-scale processes are thought to he the primary proc- 
esses for cyclogenesis. frontogenesis is a secondary process that may feed back on the 
synoptic-scale forcing to produce stronger cyclogenesis than that given by just the 
synoptic-scale forcing. The interaction between synoptic-scale vorticity and temperature 
advection and mesoscale frontogenesis in cyclones is not yet completely understood. 

1. Sub-Synoptic Scale Processes 

Recent ccclogenesis studies suggest important coupling between cyclogenesis 
and mesoscale processes associated with frontogenesis and jet streaks. L'cccllini et al. 
(19S5. 19S7) and Keyscr and Shapiro (1986) have shown that small scale upper-level 
baroclinic processes associated with upper-level jet streaks enhance the vertical circu- 
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lation and cyclogenesis. Strong jet streak divergence induces a secondary circulation 
that extends through the depth of the troposphere. This enhances the development of 
a low-level jet and cyclogenesis at the surface. This mesoscale process is considered to 
be an important mechanism for some explosive cyclones. 

Another important mechanism in some explosive cyclones is the release of latent 
heat in the mid-troposphere. The release of latent heat also influences the frontogenesis 
process. Williams et ‘al. (I9SI), using a numerical model, and Thorpe vl984). using an 
analytical quasi-geostrophic model, showed that latent heat release on the warm side of 
the front is frontogenetical aloft. The heating directly increases the temperature gradient 
at levels of large heating, and the induced divergent circulation leads to convergence 
below the heating which increases the temperature gradient in that region. The precipi- 
tation maximum occurs where the vertical motion is a maximum. Evaporative cooling 
occurs when precipitation falls through dry air below on the cold side of fronts, thus 
enhancing the cross-front temperature gradient ahead of the front. 

2. Planetary Boundary Layer Processes 

.Although baroclinic synoptic and mesoscale processes mentioned above are 
dominant, the importance of surface and boundary layer processes also has been dem- 
onstrated in recent studies. The boundary layer processes contribute to cyclogenesis and 
frontogenesis through thermodynamic effects due to the heat and moisture fluxes, and 
dynamic effects due to friction. Friction further disrupts geostrophic balance and in- 
creases the ageostrophic component of the wind through frictional wind turning. This 
results in an additional component of convergence and shear due to friction. Keyser and 
Anthes (1982) showed that convergence due to surface friction acts frontogenctically at 
the surface. Behind the frontal zone the vertical shear of the ageostrophic wind, which 
is frictionally induced, acts to diminish the static stability in the planetary boundary 
layer (PBL) and enhance the static stability at the top of the boundar}' layer. Ahead of 
tlie front, differential temperature advection by the ageostrophic component of the wind 
acts to increase the static stability in the boundary layer and decrease the static stability 
at the top of the boundary layer. The result is a compressed thermal gradient in the 
PBL. which is frontogenetical. Keyser and Anthes (1982) also showed that frictional 
convergence dri\es a vertical velocity jet above the surface front. This vertical jet is a 
component of the secondary circulation. 

Danard and Ellenton (19S0) evaluated the surface heating in several cast coast 
U. S. cyclones, and concluded that the surface heating acted to decrease the low-level 
barocliniuty by warming the cold advection regions and vice versa. This finding is 



consistent with previous studies of the eflects of surface heating in baroclinic instabilit}' 
models, and indicates a frontolytical role of the surface heat flux. However, in studies 
of the Presidents' Day storm of 1979. Bosart (1981) and Bosart and Lin (1984) indicated 
that the low-level contributions by the surface heat and moisture fluxes were important 
to cyclogenesis. They found that positive heat fluxes to the northeast of the surface low 
potentially contributed to rapid development through distribution of the cyclogenetic 
environment. Reed and Albright (1986) found similar results in their study of an eastern 
North Pacific cyclone. 

A favorable surface heat and moisture flux distribution can also enhance 
frontogenesis and cyclogenesis through the associated latent heat release above the 
boundary layer. Chen et al. (1983) found that mid-tropospheric latent heating decreased 
when the surface heat and moisture fluxes were removed. Atlas (1987) found that the 
vertical transport of air that is heated and moistened by the surfitce heat and moisture 
fluxes can result in significant mid-tropospheric heating through the release of latent 
he.it. 



B. MARINE CYCLONES 

Compared with continental and coastal cyclones, marine cyclones have a more 
complex surface heat and moisture flux distribution contributing to frontogenesis and 
cyclogenesis. due to the influence of the sea-surlace temperature (SST) distribution. 
Busingcr and Shaw (1984) suggested that for a simple uniform flow across an SST gra- 
dient. PBL convergence will result from stratification eflects on the frictional wind 
turning as the wind flows across an SST gradient. As a cyclone moves across a fixed 
SST distribution, changes in the frictional convergence and the transport of boundary 
layer air to the mid-troposphere are likely to occur. In their study of surface fluxes and 
divergence in Storm Transfer and Response Experiment (STREX) storms. Fleagle and 
Nuss (1985) showed that the distributions of the frictional convergence and vertical ve- 
locity at the top of the PBL are primarily determined by the horizontal distribution of 
the surface wind stress. In turn, the wind stress depends not only on the synoptic-scale 
forcing at the top of the PBL. but also on the PBL stratification, the surface heat flux 
distribution, and thermal wind turning in the PBL. Nuss (1989) showed that stable 
stratification behind the warm front, and unstable stratification ahead of the warm front 
in a region of upward heat fluxes, apparently resulted in enhanced frictional convergence 
ah.ead of the front. Danard (1986) found that increasing the SST by five degrees resulted 
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in greater boundary layer convergence and enhanced convective heating in the model 
cyclogenesis, 

Nuss and Anthes (1987) studied the ellcct of surface fluxes and boundary layer in- 
fluences on an idealized marine cyclone. They showed that although baroclinic insta- 
bility, which depends on the \ertical wind shear and static stabilit}', as shown originally 
by Charncy (1947) and Eady (1949), is sufl'icient for rapid cyclogenesis, relatively small 
changes in the low-level static stability, meridional temperature gradient, and surface 
heat and ntoisttirc fluxes significantly modified the development of the model cyclone. 
These low-level changes resulted entirely from varying the SST distribution. They found 
that the deepening rate of the cyclone increased by 25To for complete removal of the 
surface heat and moisture fluxes, and increased by 157o for a modified SST distribution 
in phase with the low-level temperature pattern. 

These results were further analyzed b}' Nuss (1989), in order to clarify the role of 
surface fluxes in the overall development of the cyclone. lie found that the effect of 
surface heat and moisture fluxes on the structure and dynamics of oceanic ct'clones can 
be substantial and that strictly boundary layer processes act to enhance the vertical cir- 
culation and cyclogenesis through a dynamic feedback process. Specifically, he found 
that the cfTcct of surface heat and moisture fluxes depends critically on their distribution 
rclatise to the cyclone. Tor a zonal SST distribution, patterns of surface heating 
counteract I’BL temperature advection and oppose cyclogcnesis, as found in previous 
studies, including Danard and Ellenton (1980). Tor a sinusoidal SST distribution es- 
sentially in phase with the low-level temperature distribution, the resultant pattern of 
surface heating yielded unstable stratification to the northeast and stable stratification 
south of the warm front. This acted to enhance cyclogcnesis through enhanced frictional 
con\ ergence as well as increased heat and moisture transport to the surface cyclone from 
the northeast. 

The relationshiji of enhanced frictional convergence to other frontogenesis processes 
was not examined and is the focus of this study. Although the surface heating pattern 
may be a damping effect, associated mid-tropospheric latent heat release increases the 
temperature gradient and cyclogenesis, which indirectly enhances the PBT frontogenesis 
processes. This indirect relationship between surface diabatic processes and the intensity 
of the gcostrophic frontogenesis and cyclogenesis is not well understood. Ucccllini et 
al. (19S7) have demonstrated the importance of the non-linear coupling between 
botindar}' layer processes, latent heat release, and jet stream forcing along a coastal front 
for the Presidents' Day storm. Tor open ocean cyclones, the warm front region is likeh’ 
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to be the focal point of the interaction between surface fluxes, diabatic heating, and the 
jet. as suggested by Nuss (19S9), and will be the subject of this study. Studies by 
Emanuel (19S5, 198S) and Sanders (19S6) also suggest that latent heat processes associ- 
ated with the warm front may contribute to enhanced cyclogenesis over the ocean. 

C. OB.IECTl\E 

The purpose of this research is to quantify the contributions of various forcing 
mechanisms to the secondary circulation around a warm front, and assess their impact 
on the development of a cyclone. The influence of the surface heat and moisture fluxes 
on the frictional forcing and diabatic heating, and resultant secondary circulation, is of 
particular interest since it is not as clearly understood as the larger-scale geostrophic 
deformation processes. In addition, the feedback of these diabatic processes on the 
geostrophic deformation and warm frontal intensification will be examined. The specific 
objectives arc: 

1. Diagnose the geostrophic versus diabatic processes in warm frontogenesis 
to distinguish frictional convergence from purs geostrophic deformation; 

2. Describe the feedback mechanism of surface heat fluxes on geostrophic and 
diabatic forcing by quantifying their frontogcnetic effects on the friction 
and latent heat release and their frontolytical effects on the PBL thermal 
gradient; 

3. Describe the time evolution of the forcing; 

4. Obtain an understanding of the interaction of the SS f tvith fronts. 

The Sawyer-Eliassen diagnostic equation will be used to diagnose the forcing of the 
vertical circulation around the warm front in three numerical simulations of 
cyclogcnesis. The Sawyer-Eliassen equation is reviewed in Section II. followed by a brief 
description of the model experiments in Section III. Section IV compares the results 
of the three simulations. A discussion of these results is given in Section V, followed by 
conclusions and recommendations for future research in Sections VI and \TI. 
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II. BACKGROUND 



Sawyer (1956) showed how the secondars' circulation around fronts can be calcu- 
lated on the basis of quasi-geostrophic theory, and Eliassen (1959, 1962) extended his 
work to include friction and latent heating effects, as well as the vertical shear in the 
cross-front plane, to the forcing of the secondary circulation. They developed the theor}' 
of describing temporal changes in the vertical and horizontal gradients of temperature 
and absolute momentum in the cross-front plane in terms of the induced geostrophic and 
agcostrophic deformation field. It is from this theory that the Sawwer-Eliasscn diag- 
nostic equation for the secondary circulation was derived. This equation enables diag- 
nosis of small-scale diabatic and frictional processes as well as those processes associated 
with geostrophic deformation forcing the secondary circulation around a front. Eliassen 
(1962) demonstrated that the circulation is counter-clockwise (cyclonic) around regions 
of positive forcing, and clockwise (anti-cyclonic) around regions of negative forcing. The 
circulations are concentric ellipses in the vicinity of the point source of forcing (Fig. 1). 
oriented vertically along the front and becoming vertically or horizontally distorted at 
some distance from the point source as a function of vorticity and inertial stability 
iKeyser and Shapiro 19S6). The Saw}er-Eliassen equation was later expanded by 
Hoskins and Draghici (1977) and Hoskins et al. (197S) to account for variations in the 
.along-liont plane, thus generalizing the theory of secondar}' circulation into three di- 
mension^. known as the Q-vector theory. The resultant equation is known as the 
three-dimensional (3-D) form of the Sarryer-Eliassen equation. 

A. ACEOSTROPHIC C IRCULATION AROUND A WARM FRONT 

Gidel iI97S> describes the processes of geostrophic balance and continuity of mass 
flow which lead to the secondary circulation around a warm front. The rising motion 
within the front coupled with the sinking motion on the cold air side of the front is a 
thermodynamically direct circulation. By continuity of mass, the flow is convergent at 
low levels in the upward vertical motion region and divergent at upper levels. Transverse 
flow results at lower and upper levels, completing the circulation cell in the region. 
Upper-level westerly (positive) momentum and low level easterly (negative) momentum 
are created, which contribute to increasing the vertical shear of the along-front wind 
component. The adiabatic cooling of the rising warm air and adiabatic warming of the 
sinking cold air contribute to decreasing the temperature gradient across the front. Figs. 
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Fig. 1. Streaiiiliiies aiouiul a positive point source. Streamlines in the (a) /n/;-phine 
and in the (b) jy’-jrlanc. m denotes absolute moincnluni. From Fliassen 
(1962). 



2 and 3, from Gidcl (1978), depict the initial and advected potential temperature surfaces 
and isotachs, and the induced secondary circulation around a warm Iront. In the frontal 
zone the wind field is stronger than that required for geostrophic balance, implying that 
Coriolis dominates the which sets up a transverse flow to the right (north). Im- 

mediately ahead of the front, the opposite is true, and the resultant transverse flow is to 
the left (south). I'he flow is convergent directly beneath the front, and divergent behind 
and ahead of the front. 1 he convergence near the leading edge of the frontal zone acts 
as a deformation field which enlKinccs the potential temperature gradient across the 
front. The circulation cell behind the front is thermodynamically indirect and is the in- 
direct cell discussed by Fliassen (1962) with respect to warm fronts. In the indirect cir- 
culation, the cross-front thermal gradient is enhanced by adiabatic processes and the 
along-front vertical wind shear is diminished by low-level divergence and upji .r-lcvel 
convergence. 



7 




Fig. 2. Idealized cross section through a \\arni front. l]\e initial position of the 
potential temperature surfaces is indicated hy the solid lines, and their po- 
sition after they have been slightly advected is indicated by dotted lines. 
1 he dashed lines are isotachs of the along-front wind held in geostrophi’e 
balance with the initial potential temperature field. From Gidel (1978). 
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Fig. 3. Secondary circulation around a narin front. The dotted lines are repeated 
from Fig. 2. The daslicd lines are the isotachs ol'tlie along-f'ront wind held 
after being advected from their initial position shown in Tig. 2. 1 he solid 
lines are the isotachs of the along-front wind field u which geostrophically 
balances the potential temperature field, shown as dotted lines. The heavy- 
arrows indicate horizontal and vertical components of the induced circu- 
lation. From Gidel ( 1978). 



Ihe ageostrophic flow around a positive center of forcing is counter-clockwise, 
which gives the sense of the secondary direct circulation cell if it was located ahead of 
the warm front. Similarly, a negative (clockwise) center gives the direction of flow of the 
indirect cell if the center was located behind the front. Fig. 4a is a depiction of the the- 
oretical location of the positive and negative ellipses of forcing around a warm I'ront. 
A dipole effect is set up between the two centers, and the gradient of forcing from one 
cell to the other is depicted as strongest across the vertical flow regions, assuming the 
downward and upward motions are the same magnitude. '1 heoretically, the stronger the 
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magnitude and gradient of the cells, the stronger the magnitude of How of the circulation 
cell. If the updraft region is stronger than the downdraft region, one would expect the 
gradients to be as in Fig. 4b. 

The above discussion assumes the scale of geostrophic and diabatic forcing is the 
same as that of the circulation cell. However, processes which contribute to enhance- 
ment of the circulation cell are not on the same scale. The Sa^^wer-Fliassen equation is 
a means by which the small-scale processes of friction and diabatic heating can be in- 
cluded with geostrophic deformation to diagnose the forcing of the cell. Total forcing 
would not necessarily be a single symmetric center, but rather a center distoited by 
smaller regions of positive and negative forcing acting to enhance or retard the second- 
ary circulation. For example, one would expect enhanced regions of forcing at low levels 
due to geostrophic and frictional conlluence, and at mid-levels due to latent heat release. 




Fig. 4. Ageostrophic forcing centers around a warm front. L'llipses of forcing are 
for (a) an updraft and downdraft of equal strength, and (b) a stronger up- 
draft. Note the orientation of the front and x- and y-a\es have been re- 
versed from Figs. 2 and 3. 
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B. SAWVER-ELIASSEN EQUATION 

Eliassen (1962) introduced the quantity absolute momemwn to describe the second- 
ary circulation in terms of an absolute frame of reference, the equations for which can 
be written as 



oK^u.-iy ( 1 ) 

and 



HI, = w +./:v 

The thermal wind relationships in pressure coordinates are 






''r 
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( 3 ) 



and 
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following Hoskins and Draghici (1977), and Keyser and Shapiro (1986). prognostic 
equations lor absolute momentum and thermal advection can be written as 
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where JJc,i{) = -7^-Ct “ ‘Vt'^ is the Jacobian operator, 
cc cd cd cc 
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/•_ = — — . and F,, = — — , where 1 is the stress on a horizontal surface, 

( P cp 
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By rearranging the terms, the x and y components of the Sawyer-Eliassen equation 
can be written as 



U„{'V.,) + ycpo- cy) + ( -^ ^ )«]/ = [ -2J„{„,V,) + ^ - V ^ (9) 



and 
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By taking the dot product of — A X yV6 and the vector equation made up of the x 
and y components given above, the magnitude of the forcing and circulation can be 
obtained. The right hand side (forcing) of the equation becomes 
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Tlie terms are then normalized by di\iding by | /< X y'^0 |. to obtain the magnitudes 
of the components contributing to total kinematic forcing of the warm front. The 
Sawcer-Lliassen equation is not solved; rather this form enables a diagnosis of the con- 
tributions to the forcing of tlie circulation cell b> each term. Centers of positive and 
negative forcing arc derived by this diagnosis, and represent centers of counter-clockwise 
and clockwise circulation respectively. The 3-D form is used in place of the 2-D version 
because the data set is three-dimensional and because the along-front gradients are not 
assumed to be in gcostrophic balance. 



C. DCSCRIPTIO.N OF TERMS 
1. Geostrophic Confluence 

The first term in equation (11) represents the effect that confluent gcostrophic 
deformation has in enhancing the horizontal temperature gradient (V^J') b\' compressing 
the isotherms across the front. The strongest centers of confluence (positive forcing) and 
diniuence should be as in Tig. 5. The low-level confluence and upper-level diflluence 
represent a downward transport of momentum, which acts to reduce the vertical wind 
shear, while at the same time the thermal gradient is increasing. Following Keyser and 
Shapiro (19S6). the resultant contributions to 4"^~) / and d[y-^) / dt arc 'espec- 
ti\'cly negati\e and positive, which is a tendency- towards imbalance in the thermal wind, 
a frontogcnetic process. This imbalance is counteracted b>' the secondare’ direct circu- 
lation lorccd bv the confluence, which simultaneouslv increases the vertical wind shear 
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tlirough tlie net upward transiiort of positive momentum, while reducing the temper- 
ature gradient by adiabatic processes. 

2. Geostropliic Shear 

The second term represents the elTcct that geostropliic shear deformation has in 
rotating the along-front tcniperature gradient perpendicular to the front. If the winds 
across the front increase to a maximum at the low level easterly jet on the cold air side 
(ahead) of a warm front, cyclonic shear is present. If the easterly wind ahead of the front 
is stronger than that behind the front, dilTerential cyclonic shear in the presence of 
along-front decreases in potential temperature tends to enhance the cross-front temper- 
ature gradient. Larger negative values of u, ahead of the front relative to weaker values 
c 

behind the front dilferentially advect higher values of potential temperature 

at a greater rate north of the front relative to south of the front. At the same time, the 

\ ertical wind shear (-^) tends to tip the u. isotachs into the vertical where the horizontal 
cp 

shear is cyclonic." thus reducing the magnitude ol the vertical wind shear, creating an 
imbalance. Again, the direct circulation would counteract the imbalance created by the 
shear as with the confluence. One would expect that the strongest forcing due to shear 
would be between the low-level easterly jet and the front. 




Fig. 5. Transverse ageoslrophic circulations foi- an idealized front. Circulations 

dr„ cv, 

characterized bv (a) conducnce ( < 0 , - 7 ^ = 0 ) and (b) horizontal shear 

cv^ cVj ' O' rp 

( — = 0, — > 0). Dashed lines, isotachs of u. (denoted bv U); thin solid 

vy op 

line's, isotachs of (denoted by V); thick solid lines, streamlines of 
ageostrophic circulation. Adapted from Tliassen ( 1962 ). From Keyser and 
Shapiro (1986). 
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3. Friction 



Tlie tliird term represents the convergent eflect that friction has on the 
geostrophic wind. Winds at the surface experience friction, but are maintained in part 
by the downward momentum flux from stronger winds above the surface. DiHerential 
downward momentum fluxes act like confluence to compress the temperature gradient. 
.As the cyclone intensifies the winds increase, which leads to increased friction and fric- 
tional turning of the winds toward the low. This enhances frictional convergence. If the 
winds were to increase uniformly around the low, there would be no differential mo- 
mentum flux across the front. However, the wind gradient across a frontal zone is 
greater than elswhere around the low. so the maximum differential friction is in the vi- 
cinit}' of the fronts. 

As discussed in Section I. favorable surface heat fluxes also enhance frictional 
convergence. Regions of upward heat flux into the PBL from a relatively warm surface 
can be likened to warm air advection into the region. This would cause a veering of the 
PBL wind with height, or an increased turning angle and further convergence at the 
surface. Similarly, a downward heat flux from the PBL to a cool surface would cause a 
backing of the thermal wind with height, or a reduced angle of turning and reduced 
convergence. Lor a warm front in which there are downward heat fluxes in the cold air 
sector ahead of the warm front as well as in the warm air sector behind the front, one 
can expect that for a given SST distribution, the air-sea temperature diflercnce would 
be greater in the warm air region behind the front than in the cool air region ahead of 
the front. '1 he downward heat flux would be greater behind the front than ahead, with 
the backing of the thermal u ind stronger behind the front than ahead. This creates a 
net din’erential across the front, with stronger convergence ahead of the front. Thus, the 
f rictional convergence is composed of the geostrophic component and a heat flux com- 
ponent. I ig. 6 depicts the geostrophic frictional forcing as a result of the difTerential 
downward transfer of momentum from tlie winds above the PBL. and the heat flux 
component of frictional forcing within the PBL due to surface heat fluxes. The depiction 
is for the case of stronger downward momentum fluxes behind the front than ahead, and 
for stronger downward heat fluxes at the surface behind the front. This would result in 
a reinforcement of frictional confluence. The greatest difTerential would be across the 
^^■arnl front, and the strongest frictional confluence would be located where the strongest 
differential occurs. Realistically it is possible to have stronger downward momentum 
fluxes ahead of the warm front, which would oppose the heat (lux component depicted 
abo\e. resulting in negligible frictional forcing. 
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rig. 6. Cenlers of fiidional foixing. 1 Ica \7 lines above PBL denote secondary 
circulation cells, with vertical arrows in the PBL denoting downward mo- 
mentum lluxes (solid) and downward heat (luxes (curved). Ileavw arrows 
indicate stronger downward momentum lluxes. Dashed lines depict warm 
liont. 

4. Diabatic Heating 

The fourth term represents the effect diabatic healing has on the horizontal 
temperature gradient. In the boundary layer, dilferential heat fluxes act to relax the 
temperature gradient, and are strongly dependent upon the surlace temperature distrib- 
ution, as discussed in Section 1. Surface heat and moisture lluxes in the PBL can be 
turbulently mixed out of the PBL. providing a source of heat and moisture to the mid- 
troposphere in addition to the heat and moisture advected into the region. At mid- 
levels, latent heat release within and behind the warm front, and evaporative cooling 
ahead of the front enhance the temperature gradient. One should expect centers of 
positive and negative diabatic foicing to be located as in Tig. 7. I he low level forcing 
opposes the secondary circulation at low levels and the mid level forcing enhances the 
circulation. 
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Fig. 7. Centers of diabatic foieing. Forcing in the PHI opposes the secondary 
circulation. In tire mid-troposphere it enhances the circulation. 
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111. DESCRIPTION OF DATA SET 



A. M .MERICAL MODEL 

The numerical simulations for the three case studies used in this study were run for 
4S hours using the Pennsylvania State University National Center for Atmospheric Re- 
search (PSU NCAR) mesoscale model version 4 (MM4). which is described in detail by 
Anthes et al. (19S7). The usual meteorological data plus the boundary layer fluxes were 
output even.' 3 hours. The version of the model used for this study utilized 14 vertical 
layers from the surface to 100 mb. periodic cast-west boundaries, a /I - plane variation 
of the Coriolis parameter, and the high resolution boundar\- layer model. The boundarv' 
layer occupies live of the 14 layers, with the lowest layer within the surface layer. 

DiPcrcntial friction can result from the boundary layer structure and processes in 
the model since the model vertical heat, moisture, and momentum fluxes use a vertical 
dilTu'iion coefficient that depends on the local bulk Richardson number (R,,). The 
boundary conditions are considered stable when the bulk Richardson number (Rf) is 
greater than zero and Z^jL is large. Zj, is the height of the surface layer and L is the 
Monin-Obukhov length. Tor very stable conditions (/?,^>0.2 and unrestricted), 
no turbulent fluxes occur. For stable conditions (0.0 < R^< 0.2 and ZfJL unrestricted), 
damped mechanical turbulent fluxes occur, and for marginally unstable or neutral con- 
ditions (Ri,< 0.0 and Z,JL < 1.5). both mechanical and buoyant turbulent fluxes occur. 
Unstable boundaiw conditions apply when both R^<0A) and Zf^jL> 1.5. and under 
tftese conditions the vertical fluxes are computed using a convective plume model. In- 
cluded in the fluxes is a downward entrainment flux at the top of the ntixed layer equal 
to 20'’ n of the surface flux. 

Cloud formation and precipitation are parameterized through both non-convective 
and coTweetive processes. Stable lifting that results in relative huntidities greater than 
10(U,. is treated as non-convective precipitation, and the latent heat is added directly to 
the thermodynamic equation. Unstable ascent or convection is parameterized using a 
modified Kuo (1974) scheme (Anthes 1977) which redistributes the moisture and latent 
heat \ertically through prescribed profiles, fhe criteria for convection arc total column 
moisture convergence greater than U.Ol g m - s ' and the predicted cloud height greater 
than mb. 
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DiLibmic heating used in the Sawyer-Eliasscn diagnostic equation is computed as the 
residual of the thermodynamic equation 

O = -^ = C,( ^ + I ' . r + « -g- ) + a.., 

where 

. Pr . R 

0 = 7( — )e.\p — 

‘ '“f 

The time tendency was computed as a six-hour time difTerence and introduces some error 
into the diabatic heating values. However, the diabatic heating data does correspond 
with the rainfall regions, giving at least qualitative confidence in their distribution. Be- 
cause frictional forcing is computed as a second order vertical difTcrential. forcing at the 
surface is reprcsenied by forcing at a level above the surface in the middle of the PBL. 

B. IMIIAL C ONDITIONS 

■fhe initial sea-level pressure (SLP). surface equivalent potential temperature, and 
ebserwd wind held I'or the three cases arc shown in Tig. S. The initial baroclinic insta- 
bility conditions arc the same in each case, with the only dilTerences due to the respective 
surface heat and moisture flux distributions. The initial model cyclone diflers from an 
actual atmospheric cs clone most significantly in its lack of an initial vertical circulation, 
its moisture structure being too dry compared with typical cloud patterns of actual 
cyclones, and a slight upshear tilt in the thermal wave. All but the moisture distribution 
diflerence are eliminated in the first six hours of each model simulation. The moisture 
distribution takes IS hours to resemble the typical cloud pattern of an actual cyclone, 
which also gives less certainty to the diabatic heating contributions to the model 
cyclogenesis and frontogenesis. 

The control experiment ^Case 1) uses a zonal SS T distribution (Fig. 9) and includes 
surl'ace heat and moisture fluxes throughout the simulation, ^^'lth this SS f distribution 
the mean baroclinicity in the atmosphere is reinforced, and the pattern of surface heating 
introduces temperature tendencies that counteract the initial thermal ad\ection (Nuss 
and .Anthes. 19S"). The initial surface flux distribution is characterized by upward heat 
and moisture fluxes west of the surface low and downward fluxes east of the surface low. 
The second experiment (Case 2) uses the same zonal SST distribution but insulates the 
atmosphere from the sea surface by specifying a zero heat and moisture flux condition 
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Fig. 8. Initial conditions. Sl.P in mb (solid), 0, in K (dashed), and observed wind 
vectors in m s "' lor each case. Length ol' wind vector denotes intensity, 
with the maximum length representing 50 m 5''. 

throughout the simulation. Because the only source of moisture in the model is li om the 
surface moisture flux, this case has no moisture directly transported to mid and U]-iper 
levels as in the other two cases, and relies on the cloud formation and precipitation 
parameterization applied to rising warm air cooled through adiabatic processes, 1 he 
downward momentum llux is not specifically modilied from Case 1, except indirectly 
through stratification elfects that result from no surface heating, ^\■ith no surface Iluxcs, 
the total frictional convergence is reduced to that caused by dilferential downward mo- 
mentum fluxes, and there is no thermal damping effect in the boundary layer. 
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Fig. 9. Case 1 and Case 2 SSI distribution. SS T in °C (solid), with tlic initial total 
(latent plus sensible) lieat Ilux distribution in W /?? and the cyclone track 
and position at six-hour intervals from 00 h to 48 h indicated for Case 1. 



The liiial experiment, Case 3. uses a sinusoidal SS I' distribution and includes surface 
heat and moisture fluxes as in (iase 1 (I'ig. 10). The sinusoidal SS f distribution in this 
case is in phase with the thermal structure of the cyclone, thus reinforcing the low-level 
baroclinicity instead of the large scale mean baroclinicity. Ihis distribution results in 
large upward total heat fluxes to the northeast of the surface low as well as to the west. 
In this case the magnitude of the heat flux is largest northeast of the surface low, en- 
hancing the frictional convergence northeast of the surface low compared with Case 1. 
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Fig. 10. As in Fig. 9, except for Case 3. 

C. DEVELOPMENT OF THE ( VCLONE 

'I'he dillercnces in structure and dynamics of each cyclone throughout the simu- 
lations represent tlie direct and indirect cfl'ects of the surface flux distribution unitiuc to 
each case, since the initial barodinic instability factors are identical. The intensity of 
cyclogenesis in each case can be described in part by the sea-level pressure (SEP) 
deepening rate. Fig. 1 1 is a plot of the lowest closed contour of SEP around the center 
versus time for each case. By the Sanders and Gyakum (1980) definition of explosive 
development (1 mb /r' for 24 h or more), all three cyclones develop explosively, with the 
average deepening rate 1.46 mb /;■' for Case 1, 1.89 mb /i"' for Case 2, and 1.65 mb ' 
for Case 3 (Nuss and Anthes 1987). 

In Case 1, the cyclone tracks northeast towards colder water, initially over n region 
of increasing downward heat fluxes, and then towards the region of increasing upward 
heat fluxes in the top right quadrant of Fig. 9. As the cyclone moves towards this re- 
gion, the differential in static stability across the warm front is enhanced. Incicased 
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3 arc equivalent until 30 li. 

upward heat fluxes ahead of'llic front reduce the static stability in this region compared 
with the comparatively small increase in static stability behind the front in the region 
of downward heat fluxes. Dilfciential static stability in the PHL increases frictional 
convergence, however it is small relative to the thermal damping effects of the heat flux 
distribution, which acts to slow liontogcnensis by warming the cold air sector ahead of 
the front and cooling the warm air sector behind the front. '1 his prevents the formation 
of concentrated temperature gradients at the surface and associated strong tempciature 
advection (Nuss and Anthes, 1987). 

The removal of surface heat and moisture fluxes in Case 2 results in an increased 
rate of cyclogenesis compared with Case 1. In the absence of surface fluxes, differential 
static stabilty across the front is reduced, so frictional convergence due to differential 
downward momentum fluxes is no longer enhanced. However, the frontolytical thermal 
effect of the Case 1 surface heat fluxes is removed, and indeed. Nuss and Anthes (19S7) 
showed that the low-level thermal cold and warm front advection increased in this case 
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compared with Case 1, and the temperature gradients were stronger, indicating stronger 
cold and M’arm fronts. 

In Case 3 the ct clonc tracks northeast across the SST gradient, moving slightly to- 
wards warmer water until the 24 hour point (24 h) when the cyclone starts tracking to- 
wards colder water. This track results in an initial reduction in downward fluxes east 
of the surface low until 24 hours, when the downward fluxes begin to increase as the 
cyclone changes direction. The surface flux distribution of the first 24 hours favors a 
maintenance of the baroclinicity of the cyclone. The distribution of strong upward 
fluxes north of the warm front results in reduced low level static stability and stronger 
heating of the PBL in this region, which enhances frictional convergence northeast of the 
surface low. The frictional convergence in this case is stronger than in Case 1. where the 
upward heat lluxes northeast of the surface low were not as strong. Another significant 
difference between the two cases is that the moisture flux to the northeast of the warm 
front in Case 3 is greater than 200 W A/r- compared with near zero in the control case 
(Nuss. 1089). and there is an upward moisture flux throughout most of the warm sector 
in Case 3 compared with near zero fluxes in Case 1. 

As discussed by Nuss (1989), convergence near the cyclone center and warm front 
pl-'V a significant role in the vertical circulation of the cyclone. Because of the heat flux 
distribution, stronger wind stress results to the northeast of the surface low in Case 3 
compared with Cases 1 and 2. The curl of the surface stress is proportional to the fric- 
tional PBL con\ ergcncc (Fleagle and Nuss. 1985). 5\'ith increased convergence into the 
surface low and warm front in Case 3 as a result of the stronger heating, the PBL 
ageostrophic flow toward the front and low is strongest in this case. This results in a 
stronger low-let el transt erse component of the secondart' circulation, which in turn en- 
hances the intensity of the cold conveyor belt flowing towards the cyclone ahead of the 
warm front, which transports boundary layer air that has been heated and moistened 
by the surface licat fluxes toward the cyclone center. The result is increased vertical 
transport of heat and moisture out of the botindarx’ layer, contributing to stronger 
midtroposphcric heating and greater intensification of the cx clonc than in Case 1. In 
both ('ases 1 and 3. the strongest ascent was found to be along a line extending from the 
cyclone center southeast and along the warm front. In Case 3, both the frictional con- 
vergence and latent heat release are offset by the frontolytical effect of the surface heat 
flux distribution after 24 h. It appears that the reduction in the deepening rate of the 
c\ clone can be related to the removal of the favorable conditions of the first 24 hours. 
It remains to be shown whether the stronger frictional con^•crgcnce in Case 3 is due to 
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PBL processes or stronger gcostrophic forcing that results from other indirect effects of 
the surface lluxcs. 

Tlic role of surface heat and moisture fluxes in oceanic cyclogenesis is complex, The 
frontogenctical dfects of surface heat fluxes on the PBL stratification and frictional 
convergence are opposed by the direct frontolytical effects on the thermal gradient. On 
the other hand, turbulent mixing can transport heat and moisture out of the boundary 
layer into the troposphere where it can enhance the temperature gradient through the 
release of latent heat. This feeds back to the cyclone, since the enhanced thermal gra- 
dient enhances the secondary circulation and intensity of the cyclone, which in turn en- 
hances the winds and the geostrophic frictional convergence. It can be seen that the 
advance of the cyclone into a new distribution of differential heat and moisture fluxes 
has a critical impact on the intensit} of cxclogcnesis. Case 1 allows an analysis of 
boundary laser processes in the presence of a continuously changing cross-front SST 
distribution. Case 2 allows analysis of frontogenesis and cyclogenesis in the absence of 
surface lieat and moisture fluxes, and Case 3 allows analysis of these processes in the 
presence of an initially uniform SST gradient followed by a changing gradient after 24 
hours of de\ elopment. 

D. TOT.AL FORCING AT THE SI RF ACE 

1 he intensity of the surface warm front can be represented by the total geostrophic 
and ageostrophic forcing at the surface. Figs. 12 and 13 show the total forcing at the 
surface at 24 h for Cases 1 and 2. and Fig. 14 shows the same for ('ase 3. Recalling the 
direction of llor\ around \criica!ly-oricntcd centers of positi\e and negative forcing as 
discussed in Section II. one can determine the sense of flow in the vertical plane from 
the sign of the forcing at the surface. In each case there is a region of slightly negative 
forcing ahead of the warm front, which implies clockwise circulation in the vertical plane 
and consequently flow towards the northeast in the horizontal plane of the surface. 
Siinilarh . the regions of positive forcing imply counterclockwise circulation in the verti- 
cal plane, and consequently surface flow towards the southwest. The vertical plane of 
.fow is assumed to be perpendicular to the axis of horizontal forcing, and the sign of the 
forcing at the surface is assumed to extend vertically above that region. 

The boundaries between the positive and negative regions of forcing are of particular 
importance since they represent the transition from counterclockwise to clockwise flo\^'. 
and hence arc regions of upward or downward vertical motion. For example, in Case 
3. the boundary between the negatice region to the northeast of the front and the posi- 



tivc region within the surface front is a downward vertical motion region. Similarly, the 
boundary between this positive region and the negative region south of the surface low 
is an upward vertical motion region. This combination of downward motion ahead of 
the front and upward motion behind the front supports the direct circulation and is 
frontogenetical in that the forcing enhances the temperature gradient. The strongest 
surface flow is to the southwest, perpendicular to the axis of forcing, which is close to 
the position of the warm front. 

The total forcing is strongest in Case 2 and weakest in Case 1, which supports the 
results of Nuss (1989). who found that the Case 2 warm front was the strongest. "I'et this 
case did not have the strongest vertical circulation or cold conveyor belt, nor did it have 
the most rainfall, although it had the most intense cyclone. The structure of the vertical 
circulation and forcing above the surface is not necessarily the same as that of the sur- 
face. A clearer picture can be obtained from an analysis of the cross-sections of total 
forcing for each case. Of more importance is how each component of forcing is con- 
tributing to the total forcing at all levels, and what effoct the surface fluxes have on these 
components. 
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Fig. 12. Total forcing at the surface at 24 li, Case 1. Cieopotential height in ni 
(solid), total forcing in K 100 hir'^ day~^ (dashed), and observed wind 
vectors in ni s Wind \ectors vary in length according to speed, with the 

maximum lengtli 50 m 5 '. Contour interval of forcing is 5. Heavy lines 
denote zero contours. 
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rig. 13. As in Fig. 12. c.xccpl for Case 2. Contour interval is 10. 
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IV. RESULTS 



A. O^ERME^V 

A series of surface plots of total frontogenetical forcing through the cpcione domain 
were taken at six-hour intervals to analyze the evolution of total foTv.ing. Figs. 12 
through 14 are representative of the general pattern and relative intensities of forcing 
of the cold and warm fronts observed throughout the simulations. Although there is 
positive forcing of the warm front in each case, total forcing at the surface in Cases 1 
and 3 is strongest in the cold air sector behind the cold front, whereas in Case 2 total 
forcing is strongest in the vicinity of the warm front. This suggests that processes unique 
to Cases 1 and 3 are dominating total forcing at the surface. 

B. TOTAE FORCING OF THE WARM FRONT 
1. Surface 

As seen in Figs. 12 through 14, the pattern of forcing of the warm front in Cases 
2 and 3 is elliptical, whereas a distinguishable ellipse of forcing docs not become evident 
until 30 h in Case 1. fotal forcing reaches a maximum at 30 h in Case 3. and 36 h in 
Cases 1 and 2. Fig. 1.5 is a time scries of the total forcing of the warm front, measured 
as the maximum contour value at the approximate mid point of the surface warm front. 
This is in the location of maximum forcing, and also avoids the occluded front. Total 
forcing reaches a maximum of 10 K 100 km da\-^ in Case 1, 40 K 100 /n?r‘ da\^^ in 
Case 3. and 100 K 100 knr^ day-^ in Case 2. Comparing the SLP deepening rates from 
Fig. 1 1 with the rates of growth and relative intensities of total forcing from Fig. 15 for 
each Lase. one can see that the strongest frontogenetical forcing and deepest cyclone are 
both m Case 2, and the weakest total forcing and slowest deepening rate are in Case 1. 
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2. Cross Sections 

Cross sections ol' total forcing of the warm front were taken at six-hour intervals 
for each case, oriented perpendicular to the warm front and bisecting the region of 
maximum total f orcing of the warm front at the surface for that time. The cross sections 
depict more clearly the location of total forcing in the vertical plane with respect to the 
vertical circulation field and surface warm front. By 12 hours all three cases have broad 
regions of positive total forcing in the cold air sector ahead of the warm front, and gen- 
erally negative forcing in the warm air sector behind the warm front. The counter- 
clockwise circulation induced by the positive forcing ahead of the front, and the 
clockwise circulation induced by the negative forcing behind the front, set up the direct 
circulation cell and upward vertical motion region in the warm front, as suggested in 
Section II and as indicated in figs. 16 through 18 by the vectors of ageostrophic vertical 
motion. Figs. 16 and 17 show the cross sections of total forcing at 24 hours for Cases 
1 and 2, and Fig. 18 shows the same for Case 3. l ire orientation from left to right in the 
cross-sections is northeast to southwest, with the cold air sector ahead of the front to 
the left and the warip air sector behind the front to the right. 

The initial model cyclone docs not have an initial vertical circulation field; rather 
it evolves over the course of the simulation as frontogenetical processes occur in the 
model. The presence of a distinct updraft region, as depicted by the vectors of 
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ageoslrophic flow around the warm front, is evident by the 12 hour point in Case 3, with 
a less distinct updraft region in Case I at this time. Case 2 docs not develop a distin- 
guishable updraft region until 18 hours. One can get a sense of the intensity of the cold 
conveyor belt (CCB) from the magnitude of the ageostrophic vectors circulating around 
the direct cell in each case. The strongest CCB is in Case 3, while the weakest is in Case 
2 . 
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Fig. 16. 24-h cross section of total ageostrophic forcing, ('ase 1. Forcing is in K 

lOU day-' (solid), d, in K (dashed), and ageostrophic circulation vec- 
tors. Ageostrophic circulation sectors vary in length according to speed, 
with the maximum length corresponding to 50 m 5 ’ in the horizontal and 
50 /tb in the vertical. Heavy lines denote zero contours. Contour in- 
terval of forcing is 5. 
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Fig. 17. As in Fig. 16, except for Case 2. Contour inter^•al is 5. 
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fig. 18. As in Fig. 16, except for Case 3. Contour interval is 10. 



By tlic 18 hour point Cases 1 and 3 have distinguishable ellipses of positive total 
forcing extending vertically within and ahead of the updral't region, whereas in Case 2 a 
similar region is not apparent until 24 hours. Cores of positive forcing in the boundary 
layer are first distinguishable in (.'asc 1 at 30 h, and at 24 h in Cases 2 and 3, after hich 
the boundary layer remains the level at which total forcing is the strongest. Of all the 
cases, Case 2 develops the strongest boundaiw layer total forcing, reaching a maximum 
of 150 K 100 Anr' day~'^ at 36 h. The relative maxima and times of occurrence in Cases 
1 and 3 are 35 K 100 knr' day ' at 36 h and 60 K 100 knr' day-^ at 30 h respectively. 
Fig. 19 shows the evolution of total forcing in the I'BL for each case, measured from the 
respective cross sections as the maximum positive value in the core of positive forcing 
near the updral't region, typically at about 950 mb. Although this is a crude measure 
of the frontogenesis, it gi\cs at least qualitative evidence of how the forcing c\ olves in 
each case. Case 3 has the strongest initial forcing, but is the case in which the I’BL 
forcing first declines. Between IS h and 24 h the total forcing in Cases 2 and 3 is the 
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same, after which the rate of growth begins to decline in Case 3 while the growth rate 
in Case 2 continues to increase. 




At mid levels, Case 1 develops tlic strongest total forcing, reaching 35 K 100 
knr' day ‘ at 42 h. Case 2 dcsclops a maximum of 30 K 100 km ' day-', and Case 3 
develops a maximum of 20 K loo km ' day-' at 24 h. Tig. 20 is a sununary of the total 
forcing at mid levels, measured in the direct circulation cell as the maximum contour of 
positive forcing at 600 mh. 
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C. TOTAL FORCING AT 24 1 lOURS 

I'or continuity with Nuss (19S9). the 24-hour time was chosen for detailed analysis. 
This time is representative of tlic most raj^id phase of development, and is just prior to 
formation of the occluded i'ront. By this time in the evolution, the effect of the fluxes 
will have influenced all the terms to some degree. In all three cases there is a large region 
of positive forcing ahead of the warm front, and generally negative forcing beliind the 
front. The gradient of total forcing at mid levels is strongest in the updraft region where 
the vertical vectors are the strongest. The positive forcing in the direct circulation cell 
dominates the negative forcing, suggesting that the direct cell is the dominant cell in the 
secondary circulation of the warm front. 

The patterns of forcing of the direct circulation cell in each case, as seen in Figs. 
16 through 18, arc not synunctrical ellipses as discussed in Section II, particularly in 
Cases 2 and 3 where there is a cell of strong positive forcing in the boundary laser be- 
neath the updraft region. Case 1 develops a similar region by 30 hours. I hcrc a'c dis- 
tinct differences in the forcing at each level in the three cases, and as discussed in the 
previous section, must be a result of the differences in the surface fluxes in ea; Ii case. 
As described by Nuss (1989), different surface heat and moisture flux distributions re- 
sulted in differing rates of cyclogcncsis, frontogenesis, and intensity of the cold conveyor 
belt, 'fhere is certainly evidence of a relationship between total forcing of the warm 
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front, espccialh at the surface, and the cyclone deepening rate. The important question 
is how these flux distributions, which differed primarily near the warm front, modified 
the warm frontal temperature gradient and vertical motion in each case, as indicated by 
differences in total forcing. 

D. CONTRIBUTIONS BY EACH COMPONENT OF FORCING 

1. Overvims 

The total forcing in each case is strongest in the boundar}' layer, with confluence 
and shear the dominant contributors to total forcing at this level. Confluent forcing is 
generally concentrated beneath the updraft region, with shear extending from this lo- 
cation into the direct circulation cell at low levels. Diabatic heating is generally negative 
at this level in all three cases, concentrated beneath the updraft region but extending 
ahead of the warm front as with shear, frictional forcing is a small component of total 
forcing of the warm front. The forcing due to diabatic heating is the dominant compo- 
nent of total forcing at mid and upper levels, and generally located within and ahead of 
the updraft region. Centers of mid-level confluence and shear arc located within the di- 
rect cell ahead of the warm front. Interesting departures from these general results were 
found in each case, as described below. 

2. Case 1 

The 24-h sensible and latent heat flux distributions for Case 1 are shown in Fig. 
21. with the positions of the front and cross-section added for reference. By this time 
the initial 10o9 mb low has deepened to a 9S6 mb low, with strong positive total (latent 
plus sensible) heat fluxes west of the low and downward fluxes east of the low. This 
distribution prevents the formation of concentrated surface temperature gradients and 
strong temperature advcction. since the upward fluxes arc warming the cold air region 
and the downward fluxes arc cooling the warm air region. However this effect extends 
only to the top of the boundarx' layer (Nuss and Anthes. 19S7). 
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Fig. 21. 24-li heat flux (]islril)utions for Case 1. (a) sensible and (b) latent heat 

flux distributions in \V m Dashed lines indicate downward or negative 
fluxes. Fronts added Ibr reference. From Ntiss, 1989. 



a. Surface 

The frictional convergence and diabatic heating components of forcing for 
the surface are shown in I igs. 22 and 23. Frictional forcing is negative in the region of 
upward heat fluxes northeast of the surface low, but also slightly negative in the region 
of downward heat fluxes behind the warm front. However, the frictional forcing is 
strongly positive in the cold air sector behind the cold front where the upward heat flux 
is the strongest. 'Ihe inconsistency of' the frictional forcing relative to the sign of' the 
heat fluxes suggests that the momentum fluxes in the model depend on other factors that 
are not determined simply from the fluxes. This may be particularly important when the 
boundary layer is close to ncutial stratification. 

The sense of the agcostrophic flow' induced in the region of negative foicing 
northeast of the warm front is for surface flow to the northeast, with corresponding 
downward motion at the boundary between this region and the region of positive f orcing 
closer to the front. Similarly, the induced flow in this region of positive forcing is to the 
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southwest, with downward flow along the boundary with the negative region to the 
northeast closer to the warm front. The boundar>’ between this positive forcing region 
and the region of negative frictional forcing associated with the low and warm front is 
upward. The combined How would tend to support the direct circulation cell as dis- 
cussed in the previous section. 1 he overall pattern of frictional forcing resembles the 
pattern of total heat fluxes. Over the course of the simulation the band of positive fric- 
tional forcing expands in size, although the magnitude does not substantially change. 
This suggests that as the cyclone becomes more intense and the gcostrophic winds in- 
crease. downward momentum fluxes begin to influence frictional forcing more than they 
did initially, resulting in reduced negative frictional forcing. 
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Fig. 22. 24-h surface forcing due to friction, Case I. Forcing in K 100 knr' </ny"' 

(dashed), with gcopotential height (solid) in m and observed wind vectors 
in m s'. Contour interval of forcing is I. 
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Fig. 23. 24- li surface forcing due to diabatic heating, Case 1. Contour interval is 

5. 

F’orcing due to diabatic heating at the surface is frontolytical to the warm 
front, witli the values of forcing larger than IVictional forcing. Unlike the pattern of 
frictional forcing, the pattern of this component of forcing does not resemble the surface 
heat flux pattern, because the horizontal gradient of the heating is being considered. 
What is evident from f ig. 23 is that the surface heating gradient tends to reduce the 
low-level temperature gradient. J he boundary between the negative forcing region along 
the warm front and the positive forcing region behind the front implies downward mo- 
tion along the warm front. Consequently, there is opposition of the surface warm 
frontogcncsis and secondary circulation. Ihe diabatic heating component dominates 
any frontogcnetical efl'cct of the frictional forcing, giving a net frontolytical effect of the 
boundary layer processes at this time. 
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b. Cross-Sections 

The cross-sections in Figs. 24 and 25 depict tlic vertical extent of these 
diabatic forcing terms. Frictional (breing extends up to 900 mb in tlic model, at which 
point boundary layer stresses are zero. In Fig. 24, the sense of flow between the positive 
region of frictional forcing ahead of the updraft region and the negative region beneath 
the updraft is frontogenctical. supporting the leading edge of the updraft region. 
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Fig. 24. 24-h cross-sedioii of fiictional forcing, Case I. fl, in K (dashed), 

agcostrophic vertical circulation vectors in /ib 5“', and ageostrophic forc- 
ing in K 100 knr' ciav ' (solid). Contour interval is 1. llea \7 lines denote 
zero contours. 
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Fig. 25. As in Fig. 24, except for diabalic lieating. CoiUour interval is 5 . 

fhe boundary between the negative and positi\e regions of diabatic forcing beneatli the 
updraft region in Fig. 25 is frontolytical, and dominates frictional forcing as mentioned 
above. However, this frontolytical component of forcing is dominated by positive values 
of shear and confluent geostrophic deformational forcing (not shown} in the boundary 
layer ahead of the front. At mid and upper levels diabatic heating generates 
frontogcncsis and supports forcing of the updraft region throughout its vertical extent. 
The transport of heat and moisture from the surface to the mid and upper levels allows 
for latent heat to be released within and behind the warm front, which enhances the 
temperature gradient in the warm air sector. Because the forcing is positive it also en- 
hances the secondary circulation. Geostrophic deformation is also positive at mid levels, 
although diabatic heating dominates, especially near the updraft region. The diabatic 
forcing is most susceptible to error at upper levels where the horizontal temperature 
gradient is a minimum, so a further discussion of upper-level forcing is precluded, and 
is not the focus of this study. 
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3. Case 2 



By 24 h the cyclone in this case has deepened to 9S4 mb. Although it has the 
strongest warm front, it has about the same accumulated rainfall as Case 1 for this time, 
and also has the weakest secondan.’ circulation and CCB. The patterns of confluent and 
shear geostrophic deformation are similar to Case 1 and are presented for this case for 
reference. 

a. Surface 

Tigs. 26 and 27 depict the surface forcing due to confluent and shear 
geostrophic deformation, and is representative of the typical pattern of this forcing seen 
in Case 1 after 24 h. Fig. 2S depicts the forcing due to diabatic heating at the surface. 
With the absence of surface heat and moisture fluxes in this case, frictional forcing is 
negligible, suggesting that surface heat fluxes play a significant role in driving the 
downward momentum flux. Although there are no surface heat and moisture fluxes in 
this case, the pattern ol’ forcing due to diabatic heating at the surface is similar to that 
of Case 1. with stronger values. This suggests that either the errors in the diabatic 
heatiiig calculation are relatively large or that there are significant difl'crenccs in the 
boundary layer parameteri/ation in the model. 

The sense of how the positive and negative regions of forcing enhance or 
retard the secondary circulation is the same as in Case 1. Again, the boundary between 
the band of positive forcing near the warm front in Fig. 28 and the broad negative region 
to the northeast of the front is a downward motion region. Similarly the boundary be- 
tween this positive region and the negative region near the surface low is upward. As 
in Case 1. this is irontolytical. 
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rig. 26. 24-Ii confluent gcostiopliic forcing at the surface, C ase 2. Contour inter- 

val is 5. 
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Fig. 27. As in Fig. 26, e.\icpt for forcing due to shear. Contour interval is 5. 
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Fig. 28. As in Fig. 26, except for diabatic heating. Contour interval is 5. 
b. Cl oss-Sections 

The strengtli of tlic vertical circulation is weaker in this case than in Case 
1, as indicated by the vertical velocity vectors of Tigs. 29 through 31. Conlluence and 
shear, in Tigs. 29 and 30 rcspccti\ cly, are strongest in the boundary layer beneath the 
updraft region. In Fig. 31. diabatic heating at inid levels, in the I'onn of latent heat, is 
not as strong as in Case 1. This is to be expected since the only direct source of moisture 
in the model is from the PhL. Not only is there no frictional conlluence, there is no 
moisture transported by the CCB to surface low. In addition, the weaker secondary 
circulation in this case supports the weaker CCB found by Nuss (1989). However, this 
case is dominated more strongly by geostrophic deformation at the surface than Case 
1, which supports the findings in previous studies that geostrophic deformation domi- 
nates frontogenesis and cyelogenesis, since this case has the strongest warm front and 
eventually the deepest cyclone. 
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Fig. 29. 
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24-h cross-section of confluent geostrophic forcing. Case 2. Contour in- 
terval is 5. 
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Fig. 30. As in Fig. 29, except for forcing due to shear. Contour interval is 5. 



48 



2M XX MARCH 1986 5E^ EOT CVR 




0 500E* 02 
MA>awj?r^CTO« 



Fig. 31. As in Fig. 29, except for diabatic heating. Contour interval is 5. 

4. Case 3 

The cyclone in this case has de\ eloped as rapidly as the Case 2 cyclone, and al- 
though the surface front is not as intense as in Case 2 it has the most accumulated 
rainfall at 24 h. The heat lliix distiihution at 24 h, as shown in Fig. 32, is characterized 
by strong upward heat (luxes behind the cold front, positive heat Iluxcs northeast of the 
surface low and warm front, and a pocket of downward heat fluxes near the trailing edge 
of the center of the warm front, similar to Case 1. 
a. Surface 

I'he pattern of frictional forcing in Case 3, as shown in Fig. 33, is (juitc 
difTcrent to that of Case 1, in that almost the entire region east of the surface low is 
negative. Within this negative region is a pocket of slightly positive frictional forcing 
near the warm fr'ont. Again, the frictional forcing distribution is more complex than the 
heat flux pattern, indicating the significance of other factors. Whereas in Case 1 the 
r’cgion of positive forcing was ahead of the warm front, in this case it is located behind 
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the warm front. The resultant forcing of the flow between tlie positive and negative re- 
gions does not seem to located in the same position relative to the warm front to en- 
hance the direct circulation cell and frontogenesis. As in Case 1, this region of ]'»ositivc 
frictional forcing expands in si/e, although not in magnitude. ‘I he magnitude of negative 
frictional forcing does increase slightly, although it decreases in area. This latter result 
was not seen in Case 1, and suggests that since the two cyclones are of equal intensity 
at this time, the difference is due to the surface heat flux distribution. 

1 he 24-h surface diahatic heating in this case, as seen in Fig. 34, is similar 
to the other cases, but with a more narrow band of negative forcing associated with the 
warm front. 1 he strength of the gradient from the negative to positive regions is 
strongest in this case, supporting a stronger transverse flow and vertical flow at the 
boundaries. Just as frictional forcing appears to be less frontogcnetical than in Case 1, 
the locations of the upward and downward motion regions ap]->ear to be less frontolytical 
than in the other two cases. An analysis of the cross sections yields a clearer picture of 
this relationship. 




Fig. 32. As in Fig. 21, e.xeept for Case 3. 
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rig. 33. 24-h surface forcing due to friction. Case 3. Contour interval is 2. 
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Fig. 34. 24-h surface forcing due to diabatic lieatiiig, Case 3. ('ontour interval is 

5. 



b. Cross-Sections 

Ihc cross-section of fiictional forcing in fig. 35 reveals that the region of 
positive frictional forcing is located beneath the trailing edge (warm sector side) of the 
direct circulation cell updraft region, with negative forcing directly beneath and ahead 
of the updraft region. The resultant vertical motion between these centers of forcing acts 
to oppose the vertical component of the direct circulation cell. However, because the 
negative frictional forcing ahead of the front is strongest in this case compared with Case 
1, one can imagine a stronger transverse flow towards the front and cyclone at the top 
of the PBL. Thus frictional forcing is not enhancing the updralt region but is enhancing 
the transverse component of the direct circulation cell at low levels. 

This case has stronger values of mid level latent heat than the other two 
cases at this time, as seen in Tig. 36, and indeed, this case has the most rainfall accu- 
mulation at 24 h. 'fhis suggests that the region of strong upward moisture flu.x northeast 
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of the surface front has inlluenccd the latent licat of the mid troposphere, transported 
to the front and cyclone by the inlluencc of frictional forcing. 
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fig. 35. 24-li cross-section of fiictional forcing. Case 3. Contour interval is 1. 
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Fig. 36. As ill Fig. 35, except lor diabatic heating. Contour interval is 10. 



E. TIME SERIES OF FORCINC, 

The results at 24 h are a snapshot at one time, and the feedback of the surface 
diabatic elTects on geostrophic forcing has not been separated. A time series of the 
contribution to total forcing by eacli component at low, mid, and upper levels attempts 
to separate the influence of each forcing term as they evolve throughout the simulation. 
Forcing was measured as the maximum value of each term along a horizontal line across 
the direct circulation cell in the respective cross-sections of each case. The most am- 
biguous term at low levels is the frictional forcing term. Frictional forcing was measured 
as the maximum positive value near the leading edge of the updraft region of the direct 
circulation cell, since this would enhance the updraft region in this location. '1 his de- 
piction does not take into account the influence of the negative frictional forcing ahead 
of the warm front. Theoretically, the clockwise circulation around the negative forcing 
center is acting to enhance the low-level transverse component of the circulation cell at 
the top of the PBL where the circulation around the region is towards the front and 
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cyclone. Forcing due to diabatic heating was measured as the maximum negative value 
at the base of the updraft region, since it was generally frontolytical in each case at this 
level and location. Confluence and shear were measured as the maximum positive values 
in the direct circulation cell within and ahead of the front. The 950 mb level was chosen 
since it is the level of maximum values for most of the terms at each time. At mid levels, 
the 600 mb level was chosen to represent the forcing due to confluence, shear, and latent 
heat release. .Mid-level forcing was found to be a maximum at this level, with confluence 
and shear generally centered within the direct circulation cell, and latent heating located 
within the updraft region. Because the upper-level diabatic forcing is most susceptible 
to error, the maxima were located at 300 mb. where the center of the upper-level trans- 
\ erse component is located in each case. For reference, the SLP deepening rate is also 
provided with the 95o mb time series in each case. 

I. Boundary Layer 

The time series of forcing in the boundary layer for Cases 1 and 2 is depicted in 
Fig. a", and for Case 3 in Fig. 38. In all three cases, confluent and shear gcostrophic 
doformational forcing dominate total forcing at low levels, with confluence dominating 
shear throughout the simulations except for Case 1. where shear dominates confluence. 

1 otal forcing does not peak until 30 h in Case 3 and 36 h in Cases 1 and 2. 1 his is 

consistent with the decline in the deepening rate of the cyclone in each case at these 
times. Surlace diabatic heating is frontolytical, initially peaking at the same time as 
gcostrophic forcing. This component of low level forcing is frontolytical in all three 
ca^es. and reduces total I'orcing by about 20 ”o. Unlike gcostrophic deformational 
forcing, diabatic heating at the surface docs not have a single peak and period of decline. 
Frictional forcing is not a signilicant contributor to total forcing in (iases 1 and 3, but 
it is the first component of total forcing to deselop in Case 3. In all three cases, 
confluent and shear gcostrophic deformational forcing and diabatic heating do not de- 
velop at low levels until alter 12 h in Case 3 and IS h in Cases 1 and 2. The earlier de- 
velopment of gcostrophic deformation in Case 3 is presumably due to the influence of 
the frictional forcing that was present from the beginning. 1 he period of maximum de- 
velopment of gcostrophic and diabatic forcing is between IS h and 24 h in all three cases, 
and the period of decline is between 30 h and 36 h. 
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fig. 37. Time series of loieiiig in the PBL. (a) C asc 1, and (h) Case 2. Note tlie 
scale of forcing in (.'asc 2 is twice that of Case 1. 
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CASE 3 
CONFLUENCE 




Fig. 38. As in Fig. 37, except for Case 3. 

2. Mid and Upper Le\els 

'I'hc time scries of total loicing at mid levels is depicted in Figs. 39 and 40. 1 he 
time scries Ibr upper levels in depicted in Figs. 41 and 42. At both levels, diabatic 
heating dominates total forcing, except for Case 3 where gcostrophic dcformaiional 
forcing dominates. Mid level diabatic heating in the form ol' latent heat release is 
strongest in Case 1. where it does not peak and level off until 42 h. twelve hours after 
the deepening rate has begun to decline. In Cases 2 and 3 diabatic heating reaches a 
steady state at about 24 houis. At upper levels, there is some contribution to total 
forcing by confluent and shear gcostrophic deformation. In (.'ascs 1 and 2, diabatic 
heating at this level peaks at 30 h. then declines, whereas in ('asc 3 it remains steady 
after 24 h. 
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Fig. 39. As in Fig. 37, exccjil for 600 mb. 
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F. COMPARISON OF CAS IS 

1. Boundarj Layer 

Case 2 is similar to Ca<;c 1 during tlie first 18 hours, al'tcr wliich time the most 
significant dilFcrcnccs arc the rapid rise in gcostrophic Ibrcing in Case 2, the dominance 
of shear over coniluencc in C.'asc 1. and the larger surface frontolytical elFccts ofdiabatic 
heating in Case 2. Comparing Case 3 with Case I, the most significant difFcrenccs arc 
the stronger frictional forcing in the first 12 hours of development, the rapid rise in 
gcostropliic forcing similar to Case 2 lor the first 24 hours, and a decline in the 
frontolytical elTects of diabatic heating after 30 hours. 

2. Mid and L^pper Levels 

The most significant differences between Cases 1 and 2 arc the steady rise in 
diabatic heating in Case 1 cou]Mcd with the presence of geostrophic forcing. I here is 
no significant gcostrophic forcing at mid levels in Case 2. In Case 3, latent heating is 
evident, but it not the dominant mid level forcing as it is in ('asc 1. Both conllucnt and 
shear dcformational forcing dominate diabatic heating in Case 3. At upper leveh, the 
onl}- significant dilference is the absence of any significant gcostrophic dcformational 
forcing in Case 2. 
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V. DISCUSSION 



A. 0\ERV1EW 

The basic conceptual model of frontogenesis is that geostrophic deformation forces 
strong frontogenesis at surface and upper levels, and at mid levels frontogenesis is usu- 
ally weak and forced primarily by latent heat release. Frontogenetical forcing can be 
modified by small-scale frictional forcing and diabatic heating, with their relative con- 
tributions to frontogenesis modified by surface and boundary layer processes. 

A re\iew of the SEP deepening rates from Fig. 11 shows that for the first 12 hours 
each case deepens at the same rate, indicating that the surface flux distribution peculiar 
to each case does not have an apparent elfcct on frontogenesis until 12 h. None of the 
cases has significant geostrophic deformation enhancing the forcing during the first 12 
hours, suggesting that small-scale processes of frictional forcing and diabatic heating in 
the boundary layer arc important in setting the stage for the rate of development and 
structure ol' the cyclone and warm front. Once cyclogenesis has begun, forcing due to 
geostrophic deformation builds, which then feeds back on itsell'and dominates the forc- 
ing. In later stages, friction and diabatic heating act to modify the geostrophic defor- 
mational forcing. 

B. BOUNDARY LAYER 

1. Geostrophic Deformalional Forcing 

In their two-dimensional model simulation of frontogenesis. Baldwin et al. 
(1984) found that geostropliic deformation accounted for the large-scale features of the 
secondary circulation around a cold front, but was not the dominant forcing at low lev- 
els. Bond and Fleagle (1985) found similar results in their case study of an actual 
oceanic cold front. Baldwin et al. (1984) found that shear deformation dominated con- 
fluence in their model simulation. With respect to the forcing of a warm front, this study 
found that geostrophic deformation was dominant at low levels, with confluence domi- 
nating shear m those cases with the most intense warm fronts. Mudrick (1974) found 
that confluent geostrophic deformation was favorable for warm frontogenesis. with 
shear geostrophic deformation favorable for cold frontogenesis. 

As mentioned in the last section, confluence dominates shear in Cases 2 and 3. 
and shear geostrophic dcformational forcing dominates confluence at all levels in Case 
1. Shear deformation is not as efficient for frontogene.sis as confluence and may be a 
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partial explanation for the lower intensity of the warm front and cyclone in Case 1. The 
strong frontolytical eHect of the surface fluxes in Case 1 to reduce the temperature gra- 
dient also prevents the development of a frontal pressure trough. Lack of significant 
curvature precludes the development of strong confluence. This contrasts with Cases 2 
and 3. which develop stronger temperature gradients and more of a pressure trough 
along the front. In Case 3. which is the first case to develop significant gcostrophic dc- 
formational forcing, rhe frontogenctical eflect of the friction apparenth helps to con- 
centrate the temperature gradient more rapidly in this case. This occurs in spite of 
relatively strong frontolysis due to surface heating. 

2. Frictional Forcing 

Baldwin et al. (19S4) found that frictional forcing was the dominant forcing at 
low levels, and Bond and F'lcagle (I9S5) found that frictional convergence in the 
boundary layer accounted for 90*’ o of the observed total forcing at this level. The results 
of this study suggest that frictional forcing is not nearly as strong a component of forc- 
ing of the warm front, accounting for less than fn e percent of forcing at low levels. In 
addition, tlie pattern of frictional forcing ahead of the warm front is not as consistent in 
its contribution to forcing of the warm li ont as it is to the forcing of the cold front found 
in these other studies. It is not known whether this weaker forcing by friction in the 
warm front studied here is due to model characteristics or some fundamental difl'erence 
between warm and cold fronts, d he stronger frictional forcing and more intense cold 
iront in Case 3 compared with Case 1 suggest that frictional forcing can greatly enhance 
front'''genesis. This enhancement is partially due to indirect eflects on other more 
dominant processes. 

In addition, frictional forcing of the warm front is not consistently predictable 
as to wltether it is frontogenctical. .As shown in the results, sometimes the boundary 
between the regions of positive and negative forcing was such as to oppose the forcing 
by the other terms and sometimes it seems to enhance the forcing by the other terms. 
The development of these boundaries seems to be closely related to the position of the 
cyclone and warm front with respect to the surface fluxes, but apparently depends on 
other factors as well. Surface fluxes remain strongly positive beltind the cold front, 
wherea'' the flux distribution across the warm front is confused by the pocket of down- 
ward heat fluxes near the front and the track of the cyclone relative to the SST gradient. 

There is the suggestion of a time lag between frictional forcing and its effect on 
geoctrophic deformational forcing, in particular confluent forcing, in Case 3. 1 he initial 
devclnpment of confluent forcing is preceded by 12 hours of steady, positive frictional 
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forcing ahead of the warm front, and geosirophic forcing develops 6 hours prior to 
similar forcing in the other two cases, which have little or no frictional forcing. Simi- 
larly. the decrease in the rate of growth of the confluent (and shear) forcing is preceded 
by 12 hours of increasingly negative frictional forcing in this region. The reversal in the 
frictional forcing may be influenced by the track of the cyclone relative to the SST gra- 
dient. During the first 24 hours, the cyclone tracks towards the region of strong positive 
heat flu.xes and the frictional forcing is frontogenetical. When the esdone begins to 
move into colder water and away from the influence of the positive heat fluxes, the 
frictional forcing becomes frontolytical. 

3. Diabatie Heating 

Baldwin ct al. (1984). using a thermally insulated boundary layer in their model, 
found that boundary layer heat fluxes had a minimal effect on the synoptic-scale sec- 
ondary circulation. Bond and Fleagle (1985) included boundary layer heat fluxes in their 
study, yet lound similar results, although boundary layer turbulent heat fluxes were 
found to be Irontolytical in the boundary layer. Nuss and Anthes (1987) found that 
boundary las'er heat fluxes influenced the structure only below 850 mb, suggesting that 
surface heat flux effects are confined to the boundary layer, due to the influence of the 
inversion layer above the PBL. ‘fhis study found similar results in that forcing due to 
diabatie heating is Irontolytical and occurs at low levels only. Its magnitude is compa- 
rable to that of shear deformational forcing in Case 2 and confluent forcing in Case I, 
and in Case 3 it is about half that of shear. Of interest is the presence of negative forcing 
in the PBI. in Case 2 despite the absence of surface heat and moisture fluxes. This sug- 
gests that in the absence of diflerential surface fluxes, there is surface cooling perceived 
by the model. This may be the w ay the model handles convective versus non-con vective 
precipitation, or it may be a result of the computation of this term of forcing due to 
advection. 

C. .MID AND UPPER LE^ ELS 

1. Geostropliic Deformation 

As with studies of cold front forcing, geostropliic deformational forcing of the 
warm front at mid levels is dominant within the direct circulation cell ahead of the front, 
and at upper le\els it enhances the upper-level transverse component of the secondary 
circulation. However, the overall lalues are typically about 25'V of the low-level 
geostropliic forcing, 'fhe weakest mid- and upper-level geostropliic forcing is in Case 2. 
suggesting that mid- and upper-level geostrophic forcing is not as critical to the dcvel- 
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opment of the cyclone as it is at low levels, since this case has the deepest cyclone and 
strongest geostrophic forcing at the surface. 

2. Diabatic Heating 

Diabatic heating in the form of latent heat release was found to be the dominant 
forcing mechanism for the vertical component of the secondary circulation, with positive 
centers located within and slightly ahead of the updraft region and warm front. All three 
cases develop some inid-level diabatic forcing due to latent heating in iic model. The 
strongest mid-level latent heating developed in Case 1 and it continued to increase until 
42 h. This case was the weakest cyclone, which indicates that latent heating and a strong 
vertical circulation in the warm front are not the only factors for strong cyclogenesis. 
Case 3 developes a relatively steady-state of latent heating by 24 h. and the strongest 
vertical motion, however the diabatic forcing was less than Case 1. T his is presumably 
a feature of the cumulus parameterization in the model. Case 2 has the deepest cyclone 
by 24 h but produced the weakest nud-level diabatic forcing. I'his lack of latent heating 
is probably a reflection of the lack of moisture input from the surface in this case. These 
results suggest that details of the cumulus parameterization and moist processes in the 
model strongly alter the intensity of the mid-level diabatic forcing of frontogenesis and 
the secondary circulation. Since the diabatic forcing is dominant at mid and upper lev- 
els. the circulation at these levels will be sensitive to model errors in precipitation proc- 
esses. 
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VI. CONCLUSIONS 



'I he results of the model simulations conlirm the results of previous observational 
and numerical studies of front ogenesis. They also suggest important interactions be- 
tween the surface heat and moisture flux distributions, the forcing of the warm front 
secondary circulation, and cyclogenesis. Geostrophic deformational forcing was domi- 
nant. but evolved dilTerently depending upon surface friction and diabatic heating influ- 
ences. 

The presence of frictional forcing enhances the low-level transverse component of 
the secondary circulation, which in turn enhances the cold conveyor belt. With strong 
positive heat fluxes in a favorable location, as in Case 3, frictional forcing appears to 
enhance development of geostrophic deformational forcing at an earlier stage of 
cyclogenesis, which also enhances the secondaiw circulation. As the cyclone and 
geostrophic deformational forcing intensify, the surface forcing due to friction acts only 
to modify geostrophic forcing. The development of positive frictional forcing in the 
latter stages of cyclogcnesis. as in Case 1. had little effect on the deepening rate, except 
perhaps to maintain the geostrophic forcing at such a level as to yield the final 4 mb 
drop before cyclogcnesis finally ceased. 

The eflcct of surface heat and moisture fluxes and their distribution on forcing of the 
secondary circulation and cyclogcnesis is more complex, and the results less definitive 
than with frictional forcing. Surface forcing due to diabatic heating is consistently 
frontoh tical in each case. The fact that it is a strong component of total forcing in the 
absence of surface heat and moistrure fluxes, as in Case 2. suggests that either the error 
is large or there are significant boundary la}er differences in the model. Early surface 
relaxation of the frontal temperature gradient as in Case 1 appears to reduce the poten- 
tial deepening rate of the cyclone, and dominates frictional forcing. 
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VII. RECOMMENDATIONS 



The results suggest areas for further research, including: 

Verify the magnitudes of the various forcing terms in an actual case study. Because 
the model is drier than real oceanic cyclones, forcing due to diabatic heating effects may 
be different. 

Verify the effect of surface fluxes and boundary layer structure on the frictional 
forcing of the M'arm front. 

Better understand the role of various model parameterizations involved in deter- 
mining convective versus non-convective release of latent heat. This may help to explain 
the presence of negative surface diabatic forcing in the absence of surface heat and 
moisture fluxes. 
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